
###hand entered from Israeli Central Bureau of Statistics
years <- c(2009, 2014, 2019, 2024, 2029, 2034, 2039, 2044, 2049, 2054, 2059)
all <- c(1, 8, 17, 26, 36, 45, 56, 67, 80, 92, 107)
non_uo_jews <- c(1,6,11,16,20,25,29,34,39,44,49)
uo <- c(1,21,47,77,112,152,198,250,310,377,453)
arabs <- c(1,12,24,37,51,65,79,92,106,120,135)

projections <- as.data.frame(cbind(years, all, non_uo_jews,uo, arabs))

pdf('output/Figure_1.pdf')
  par(bty = 'n')
  par(mar=c(4.5, 4.5, 0.2, .2))
  plot(projections$years, projections$uo, 
       pch="", xlab="Years", ylab="Percent Expected Change",
       cex.axis = 1.5,
       cex.lab = 2,
      )
  legend("topleft", c("All","UO", "STR", "Arabs"), lty=1,
         pch=c(1, 3, 15, 24), bty = "n", cex=1.25,
      )
  lines(projections$years, projections$all, type="b", pch=1, lty=1, lwd=2,
        cex = 2.5)
  lines(projections$years, projections$uo,  type="b", pch=3, lty=1, lwd=2,
        cex = 2.5)
  lines(projections$years, projections$non_uo_jews, type="b", pch=15, lty=1, lwd=2,
        cex = 2.5)
  lines(projections$years, projections$arabs, type="b", pch=24, lty=1, lwd=2,
        cex = 2.5)
  dev.off()


